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ABSTRACT 

A  two-level  graphical -numerical  prediction  model  is  extended  so  as 
to  incorporate  some  orographical  and  non-adiabatic  warming  or  cooling 
effects  in  a  manner  such  that  almost  no  extra  work  is  added. 

The  model  is  applied  to  several  cases  where  non-adiabatic  heat 
exchanges  are  pronounced.   The  resulting  prognoses  show  a  material  improve- 
ment over  a  similar  model  which  omits  the  non-adiabatic  influence. 
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1.   Introduction 

It  is  a  well -kno-n  fact  that  during  the  winter  season  coLd  cP  air 
flowing  out  from  continent  over  the  warm  ocean  surface  will  acquire  heat 
energy  from  below  so  that  the  lower  layers  of  the  atmosphere  are  warmed. 
As  Craddock  [l]  measured  in  the  case  when  the  arctic  air  flowed  out  from 
Iceland  to  the  British  Isles,  the  heating  rate  may  be  as  large  as  65  Ay 
per  hour  and  in  the  course  of  a  day  the  1000-500  mb  layer  thickness  could 
increase  approximately  800  feet.   Hence,  in  such  cases,  an  appreciable 
error  will  be  produced  in  the  results  of  those  numerical  prediction  models 
which  assume  that  the  air  motion  is  adiabatic.   Several  techniques  have 
been  introduced  to  eliminate  this  non-adiabatic  error;  Reed  [5]  has 
incorporated  Burke's  semi -empirical  treatment  for  the  air  mass  transforma- 
tion cP-mP  into  a  graphical  model  of  numerical  prognosis.   His  results  show 
considerable  improvement  over  a  model  omitting  the  non-adiabatic  heating. 

In  the  present  investigation,  the  non-adiabatic  heating  (or  cooling) 
is  approximated  in  a  rather  crude  but  extremely  simple  manner.   As  a 
result  there  is  no  additional  work  when  non-adiabatic  heating  and  the 
orographic  influences  [3]  are  included. 


2.   Development  o£  the  prediction  equation 

The  basic  vorticity  equation  in  the  (x,y,p,t)  coordinates  takes  the 
form 


where  -— —  is  the  individual  rate-of -change  operator  following  the  isobaric 
projection  of  the  three-dimensional  motion, 

y  the  relative  vorticity, 

/     coriolis  force, 
^    del  operator. 
Both  j  and ^7  refer  to  values  measured  on  an  isobaric  surface. 
In  Eq.  (1)  the  terms  involving  the  vortex  tube  and  vertical  advection 
of  vorticity  have  been  neglected,  and  also  the  term  — fv-^y  [6]. 
Expanding  the  left  side  of  Eq.  (1), 

£(J+f)  +  V-v(f-*j)=-fv.sy  (2) 

By  equation  of  continuity, 

v-  v  =  —  ~^J- 

Eq .  (2)  becomes 

where 

Assume  now  that  to   has  a  vertical  parabolic  distribution  plus  the 
terrain- induced  vertical  motion. 


><X^.M  )  =  'ft  4#tf  VH  +C0  [,-lkte  1 


Hereto/,    0    ,  \y   are  the  surface  pressure,  density  and  wind 
velocity;  H  is  the  height  of  the  terrain;  A*  is  an  upper -pressure  level 

*  A 

(to  be  taken  as  500  tnb)  where  vertical  motion  is  strongest;  and  -Jf  is  a 
constant  of  proportionality. 

The  term  on  the  right  side  of  Eq0  (4)  involving  V/'v7H  represents  the 
terrain  induced  vertical  motion.   Here  the  surface  vertical  velocity  ^jJT  is 
assumed  to  be  proportional  to  the  horizontal  wind  and  the  gradient  of 
terrain  height  H  at  the  surface;  i.e.  &£%=  ^  \V  .  ^7  /-/    .   Moreover, 
since  (jO   — --/^  &lT"  »  ic  follows  that 


0do  =  -  gfK.X'VH 

Finally,  the  terrain-generated  vertical  velocity  is  assumed  to 
decrease  linearly  with  pressure,  and  thus  we  arrive  at  the  first  terra  const! 
tuting  :0in  Eq.  (4). 

The  second  term  of  Eq .  (4)  represents  the  commonly  used  parabolic  dis- 
tribution of  vertical  velocity  assumed  to  be  typical  of  large-scale  press- 
ure systems  [4].   Where  terrain  effects  are  negligible,  this  term  corre- 
sponds to  low-level  "isobaric  '  velocity  divergence  and  high  level  con- 
vergence (or  vice  versa),  together  with  the  level  of  non-divergence  at  p~ 

Differentiating  Eq .  (4)  with  respect  to  ^  ,  substituting  it  into 
Eq.  (3)  and  evaluating  the  result  at  A-'  gives 


^  +  V'V(l+i)  =  -  ±Jklyo.VH-*h±  (5) 

dt  £  ft 

3  * 


The  next  step  consists  of  utilizing  the  first  law  of  thermodynamics 
to  provide  a  second  equation  involving  O0m   •   The  non-adiabatic  heat 
addition  to  the  surface-500  mb  layer  will  be  assumed  to  be  proportional 
to  the  rate  at  which  the  surface  air  parcels  are  exposed  to  a  changing 
temperature  of  the  underlying  surface. 

Expanding  the  individual  change  of  the  potential  temperature,  £7   . 
and  dividing  through  by  Q      , 


The  subscript  p  at  the  ^7  operator  will  be  omitted  henceforth, 
It  can  be  shown  from  Poisson's  equation 

looo   \S" 
T 


(6) 


*  -  t(^) 


that  with  pressure  held  constant 


where  ^    is  the  specific  volume  and  R,  the  universal  gas  constant.   Sub- 
stituting these  identities  into  Eq.  (6)  yields 


V 


From  the  hydrostatic  equation, 
substitute  into  Eq .  (7)   after  multiplying  through  by  ^  ,  giving 
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It  follows  that 


£(§§) = -  V-V(%)  -  ffco  -A4$>  a>  A  &  > 


ft  V      VK-ZpJ        O^J-p-^ 


26  d  f 


(8) 


Equation  (8)  states  that  the  local  change  of  the  rate  of  change  of  the 
pressure  height  with  respect  to  pressure  is  the  result  of:   a)  advection, 
b)  vertical  motion,  and  c)  the  non-adiabat ic  heating  process  represented 
by  the  last  term. 

At  this  stage,  the  non-adiabat ic  addition  of  heat  will  be  assumed 
proportional  to  the  expression  ^  '  V  7(   where  \J/  is  the  wind  at  the 
surface  and  T   is  the  surface  temperature  of  underlying  ocean  surface. 
Hence  Eq .  (8)  may  then  be  expressed  in  the  form 

JL  (ti  \ 


where  ^  is  the  constant  of  proportionality.   Strictly  speaking,  "A\  is  not 

a  constant  but  may  possibly  vary  with  synoptic  situations  and  also  vsr) 

slowly  with  time.   However,  "rf    will  be  assumed  constant  for  the  present 

investigation. 

Now  substitute  Eq ,  (A)  into  Eq .  (9)  and  make  geoatrophic  approximation  in 

the  advection  term. 

By  integrating  between  /*•  and  /t-  ,  assuming  CS    to  be  a  constant,  we  obtain 

to  Is 


p 


Here  h  is  the  thickness  be 


tween  /i"and  y^  .  t 


In  the  last  term  of  the  preceding  equation,  \y  and  T   are  independent  of 
pressure,  hence 


g^M  +&co„-  W  iltvH  f  £ 


y,-vTs         co) 


where  if,  is  a  new  constant. 
Solving  Eq.  (10)  for  M  , 


""^  (I? +^-^  *  ^^^^H-^-V^ 


and  substituting  it  into  Eq .  (5)  gives 


-  ^Wfi+K-rt+frfcrA**4-  £  K-vTi 


Rearranging, 


<J£  *£• 


^)---v<w^+^_^      ,„, 


Nov,  we  make  the  geostrophic  approximation  for  J     in  Eq .  (11) 

o 


tfW2'"^   *'"*>/* 


$ 


According  to  Fjortoft,  the  horizontal  variation  of  -r-  may  be  neglected 
in  comparison  to  (  2),"^)  •   With  multiplication  by  — 7-3. —  ,  we  obtain 

^■vfr-i+jM&t * ^h  +igC.A.  ^M"r]  <") 


The  coriolis  force  term  on  the  right  side  of  Eq .  (12)  has  been 
manipulated  as  follows: 


Cti-3  ^ 


no,  i.t  f    ja        j-: —  *>?  -  Or 

F.Jl£££lLN 

f  ^  .—  ■  ■■■,■       ^^,— *«■ 

Therefore  Eq.  (12)  becomes 


'2. 


Jt  [z°-Z°  +C*  +G,-t  -F+NT  ] 


=  -%'V\^0~:ZB  +  Ci  +  (k  rF+A/Tj  en) 


10  v  —      l  "*i         >       I  I     I  »     /  <■      I 


« 


The  function/VTvaries  very  slowly  with  time  and  may  be  assumed 
to  be  a  constant  for  the  normal  forecast  period.   The  functions  G  and 
F  depend  on  geographical  position  only  and  are  invariant  with  time, 
hence  they  may  be  inserted  on  the  left  side  of  Eq .  (13).   Equation  (13) 
shows  that  the  quantity  (  2?p-£r   -f*  di\    "t  Qf    i    h    +  MTe   )  is  C( 


:on 


servative  with  respect  to  wind  field  at  p**  .      By  using  the  geostrophic 
approximation  for^/  ,  i.e. 


Eq.  (13)  may  be  written  in  Jacobian  form  as 


(14) 


In  order  to  arrive  at  a  more  conservative  field  for  advecting  the 
conservative  property,  the  Jacobian  identity,  J(x,y)  «  J(x+y,y),  may  be 
applied  to  Eq .  (14),  giving 

--fJfc+c4+G  +  F+rts^-%+cf+$+r+Aft]  (is) 

The  use  of  this  form  implies  that  the  [ ZQ±  C4\    ■+  Gr     +  F  +  V7^  jf  field 
is  expected  to  be  more  conservative  than  the  Z#   field  for  advection  pur- 
poses.  This  is  not  necessarily  always  the  case  since  the  thickness  field 
appears  in  the  advecting  quantity.   Nevertheless,  in  the  examples  tested, 
Eq .  (15)  gave  good  results. 
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To  depict  the  local  change  of  [_  S?e-  2^  +  Cf\  4-  Gf  ~t  p  -fN~T<  J  we 
shall  advect  this  function  for  a  short  cine  interval,  say  twelve  hours 

or  twenty-four  hours.   Denote  this  change  by  ^  .   Since  G,  P  and  NT  are 

s 

not  functions  of  time, 


=  A[M.-%+Ct]  =  -ft  (i*) 

Moreover,      -f{  ~    £    -  7 


wher 


Taking  the  space  average  over  Eq.  (17)  and  substituting  back  to  Eq.  (17), 
yields 


***=*  +  lkEB+lkAZ*] 


I       •=: 


Averaging  n   ti«es,  

v 

^=e^5^f +dr*l 


However,  this  series  converges  rapidly;  and,  as  shown  by  FJortoft  [2], 
it  can  be  approximated  as 

A%>~  6  -r  ~ ;  &  d9) 


In  the  cases  tested  in  this  investigation,  it  was  found  that  the 
second  term  on  the  right  side  of  Eq .  (19)  can  be  neglected  except  when 
a  aarked  change  is  indicated.   In  the  routine  work  where  a  time  saving  is 
of  great  importance,  /4zL  -~  £   can  ^e  used  without  much  sacrifice  in 
accuracy.   In  utilizing  the  high  speed  electronic  computer,  Eq .  (16) 
could  be  solved  to  any  desired  accuracy. 
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3.   Procedure 

In  order  to  provide  a  comparison  to  Reed's  results  [5],  the  same 
case  has  been  chosen  for  prognostication.   Tjo  other  series  were  also 
chosen  for  testing  purposes.   One  series  has  the  storm  located  in  the 
Gulf  of  Alaska.   Each  series  contains  three  consecutive  periods  with  a 
time  interval  of  twelve  hours.   Except  in  the  same  case  as  Reed's,  where 
the  prognosis  without  heating  effect  was  already  available,  two  prog- 
noses for  each  case  have  been  obtained,  one  omitting  and  the  other  includ- 
ing the  heating  effect. 

By  omitting  the  heating  effect,  Eq .  (15)  becomes 

The  following  steps  are  carried  out  in  each  prognostication: 

(1).   A  space -mean  chart  for  Z   is  graphically  constructed  using  a 

grid  distance  of  600  km. 

(2).   With  C  •  0.5,  the  current  thickness  chart  is  multiplied  by 

1/2  (this  is  achieved  by  simply  tracing  every  other  line  on  the  thickness 

chart)  and  added  to  Z  . 

o 

(3).   Next  the  result  of  step  (2)  is  added  to  the  (G  +  F)  function 
provided  by  Haltiner  and  Hesse  [3]. 

(4).   With  N  *~0  /  .  l00.J-  over  the  ocean  and  N  »  -0.05  over  the 

land  area  of  Alaska  (see  Section  4),  a  (NT  )  function  chart  was  constructed 
based  on  the  monthly  mean  ocean  surface  isotherms  published  by  the  Hydro- 
graphic  Office  and  the  current  reported  temperatures  over  the  land. 

(Fig.  19).   This  chart  of   f(NT  )  vas  added  to  the  result  of  (3).   Note 

s 

this  step  is  unnecessary  if  non-adiabatic  heating  is  omitted.   Also  it  is 

11 


suggested,  since  monthly  mean  surface  temperature  isotherms  over  the 
ocean  are  used,  a  chart  of  (G  +  F  4  NT   ocean)  function  may  be  prepared 
for  each  winter  month.   This  will  combine  step  (4)  into  step  (3)  if  the 

(NT  )  function  over  land  is  also  approximated  by  the  monthly  mean  iso- 

s 

therms  over  the  land.   Up  to  this  step,  a  chart  of  function  of 

(Z   +Ch+G+F+NT)  has  been  obtained  and  will  be  referred  to  as 
o  s 

steering  chart  henceforth. 

(5).   The  Z   chart  is  subtracted  from  steering  chart  to  obtain  the 
o 

quantity  to  be  advected. 

(6).   Next  advect  the  quantity  obtained  in  step  (5)  with  the  steering 
chart  for  twelve-hour  intervals. 

(7).   Subtract  the  result  of  step  (6)  from  the  result  of  step  (5)  to 
get  A. 

(8).   Add  \I2AZ;.    to  A  and  multiply  the  result  by  2/3  to  obtain  B. 

(9).   Average  B  to  get  B  and  multiply  by  4/3. 

(10).   Add  the  result  of  (9)  to  (8),  to  getAZ  . 

o 

(11).   Add  A  Z   to  the  initial  1000-mb  chart  to  give  the  prognostic 
o 

1000-mb  chart  of  twelve  hours  later. 
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4.   Discussion 

In  this  investigation,  three  constants  are  of  major  importance, 
namely;  C,  kj,  and  N.   Ol  these  three  constants,  k(  has  been  discussed  in 
the  Haltiner  and  Hesse  paper  [3],  and  only  C  and  N  will  be  discussed  here. 

(1).   The  constant  C. 

Several  values  of  C  have  been  used  by  various  investigators. 
In  the  present  investigation,  the  value  of  C  *  0.5  fits  best  in  the  cases 
tested.   Values  lower  than  this  will  give  an  advection  field  too  weak, 
while  higher  values  give  an  advecting  field  too  strong. 

(2) .   The  constant  N. 

In  this  investigation  the  primary  interest  is  in  the  area  of 

outflow  of  cP  air  from  inland  Alaska  to  the  sea.   Over  the  ocean  the  value 

of  constant  N  was  found  empirically  to  be  -0.1  -—gp. —      In  the  non- 

F 

adlabatic  modification  of  an  air  mass,  the  processes  which  actually  take 
place  are  not  simple.   The  degree  and  rate  of  the  modification  are  depend- 
ent upon  many  factors  such  as:   a)  the  temperature  difference  between  the 
cP  air  mass  and  sea  surface;  b)  the  roughness  of  the  sea  surface;  c)  the 
wind  speed  over  the  sea  surface;  d)  the  length  of  time  the  air  mass  has 
been  over  the  ocean;  e)  the  stability  of  the  air  mass;  f)  the  moisture  con- 
tent of  the  air  mass;  etc.   The  inter-relations  among  these  factors  are  not 
well  known  quantitatively;  consequently,  it  is  difficult  to  determine  in  a 
specific  case  the  relative  importance  of  the  various  factors.   For  example, 
air  moving  with  a  higher  speed  is  not  necessarily  being  modified  pro- 
portionately more  than  air  moving  at  a  much  lower  speed.   The  reason  for 
this  is  that  the  modification  is  not  only  a  function  of  the  distance  over 
which  the  air  has  traveled  from  the  coast,  but  also  a  function  of  the  time 

r 

the  air  mass  has  been  over  the  ocean.   The  expression  Ts  ^ '  VJf   implies 

*  °  ■> 

that  when  the  value  of  V  is  doubled,  the  amount  of  modification  will 

o 
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double,  ^hich  Is  not  necessarily  true.   Thus,  the  determination  of  the 
value  of  N  is  a  matter  of  compromise  of  the  many  different  influences. 
For  the  relatively  small  sample  tested,  the  value  of  N  *  -0.1  appears  to 
be  the  best  choice  thus  far.   However,  application  of  this  value  shows 
considerable  improvement  over  the  other  models  which  treat  the  non- 
adiabatic  heating  effect  in  a  different  way. 

In  the  course  of  this  investigation,  it  appeared  that  the  advecting 
field  tends  to  move  the  high  pressure  cell  down  stream  over  inland 
Alaska,  with  the  result  that  the  prognostic  pressure  height  was  too  high 
over  the  southwestern  side  of  the  High  and  too  low  over  the  southeastern 
side  of  the  High.   For  example,  in  the  case  1500  GMT  January  4,  1956  to 
0300  GMT  January  5,  the  high  was  built  up  on  the  east  side.   In  this  situa- 
tion, an  attempt  was  made  to  apply  the  non-adiabatic  effect  over  inland 

Alaska  by  putting  N  =  -0.05,  obtaining  a  f(NT  )  based  on  the  isotherms  on 

s 

the  current  surface  map.   Considerable  improvement  resulted.   Such  an 
approach  provides  a  continuous  temperature  field  from  the  inland  out  to  the 
sea.   This  value  of  N  was  also  applicable  to  the  other  inland  cases.   The 
combined  functions  (NTg)  for  the  land  and  the  sea  are  shown  in  Figs.  19, 
20,  21,  22.   From  these  figures,  it  seems  feasible  to  use  the  monthly  mean 
temperature  over  Alaska  land  area  as  the  basis  of  f(NTa)  over  land,  so  that 
step  (4)  can  be  entirely  incorporated  into  step  (3)  by  preparing  a  single 
chart  of/(G  +  F  +  NT  )  for  each  winter  month.   Thus  a  considerable  amount 
of  time  can  be  saved  in  daily  routine  work.   Since  the  non-adiabatic  cool- 
ing process  over  the  land,  probably  mostly  due  to  radiational  cooling,  is 
beyond  the  scope  of  this  investigation,  no  further  discussion  will  be  given. 

In  addition  to  the  empirical  determination  of  constants,  another 
aspect  at  the  advecting  stage  in  the  procedure  warrants  some  discussion. 
As  the  air  mass  moves  southward  over  an  increasing  temperature  field,  the 

14 


air  mass  will  be  warmed  up  from  below.   An  assumption  was  made  that  the 

air  parcel  directly  in  contact  with  the  sea  surface  will  immediately 

acquire  the  same  temperature  as  that  of  the  sea  surface.   In  this  case, 

the  air  mass  tends  to  be  unstable.   This  instability  will  enhance  the 

upward  heat  transport,  and  so  increase  the  modification.   On  the  contrary, 

when  the  air  mass  moves  northward  over  a  decreasing  temperature  field, 

by  the  same  assumption  it  tends  to  be  stabilized.   In  the  first  case,  the 

thickness  of  the  lower  layer  will  increase  and  in  the  latter  case  the 

thickness  will  decrease.   However  the  magnitude  in  the  first  case  of 

increasing  thickness  is  normally  much  larger  than  the  magnitude  in  the 

latter  case  of  decreasing  thickness.   Because  of  this  consideration,  an 

attempt  was  made  to  advect  less  distance  in  cold  advection  and  to  advect 

the  full  distance  in  warm  advection.   However,  this  did  not  improve  the 

results,  since  the  property  to  be  advected  is  a  combination  of  Z   Z  , 

Ch,  G,  F,  and  NT  ,  and  not  a  function  of  NT   alone.   As  a  matter  of  fact, 
s  s 

since  the  value  of  N  has  been  determined  empirically,  it  represents  a 
compromise  of  all  the  various  influences. 
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5.   Results 

The  results  of  the  tests  are  displayed  in  Table  1  which  shows  the  com- 
parison of  the  series,  with  and  without  non-adiabatic  warming  or  cooling. 
The  statistical  parameters  computed  were  based  on  a  seventy-point  grid  of 
which  forty-seven  points  were  over  the  sea.   These  grid  points  are  the 
intersections  of  the  meridians  and  latitude  parallels  which  are  5   apart 
extending  from  longitude  165  W  to  longitude  120  W  and  from  latitude  65  N 
down  to  latitude  30  N.   Due  to  the  convergence  of  meridians,  the  points 
are  denser  over  the  northern  part  of  the  area.   Since  the  area  of  denser 
points  is  inland  and  coastal  Alaska,  which  is  a  more  critical  area,  the 
values  in  Table  I  are  more  conservative  and  do  not  exaggerate  the  improve- 
ment due  to  the  inclusion  of  the  non-adiabatic  heating. 
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Figs.  1  to  16  indicate  the  results  of  four  cases.   These  are: 
Case  1)  1500  GMT  January  4,  1956  to  0300  GMT  January  5,  1956;  Cases  2, 
3,  4)  0000  GMT  January  13,  1959  to  1200  GMT  January  14,  1959,  consisting 
of  three  consecutive  prognoses.   The  other  three  cases  from  0000  GMT 
January  13,  1959  to  1200  GMT  January  14,  1959  are  not  shown.   It  has  also 
been  found  that  the  addition  of  term  NTg  gives  no  improvement  when  a 
cyclone  comes  from  the  central  North  Pacific  area;  however,  neither  is 
there  a  decrease  in  accuracy  when  f(NT  )  is  included. 

It  can  be  seen  from  the  figures  that  the  greatest  improvement  in  the 
prognostic  pressure  heights  ie  found  over  the  area  on  the  west  of  the  storm 
which  is  the  normal  location  of  the  cP  air.   This  was  true  in  all  the  cases 
tested. 

On  the  basis  of  the  sample  tested  it  may  be  concluded  that  this  model 
shows  very  good  promise.   However,  further  testing  would  be  desirable  to 
determine  the  best  values  of  the  various  constants.   The  non-adiabat ic 
heating  term  can  readily  be  incorporated  into  the  programming  of  baroclinic 
models  for  the  electronic  computer  on  which  further  testing  would  be 
greatly  facilitated.   There  is   obviously  a  considerable  saving  in  time 
compared  to  the  Reed  method  when  this  model  is  applied  manually  in  the 
graphical  technique. 
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